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We study entanglement distribution in quantum complex networks where nodes are connected by 
bipartite entangled states. These networks are characterized by a complex structure, which dramat- 
ically affects how information is transmitted through them. For pure quantum state links, quantum 
networks exhibit a remarkable feature absent in classical networks: it is possible to effectively rewire 
the network by performing local operations on the nodes. We propose a family of such quantum 
operations that decrease the entanglement percolation threshold of the network and increase the 
size of the giant connected component. We provide analytic results for complex networks with 
arbitrary (uncorrelated) degree distribution. These results are in good agreement with numerical 
simulations, which also show enhancement in correlated and real world networks. The proposed 
quantum preprocessing strategies are not robust in the presence of noise. However, even when the 
links consist of (noisy) mixed state links, one can send quantum information through a connecting 
path with a fidelity that decreases with the path length. In this noisy scenario, complex networks 
offer a clear advantage over regular lattices, namely the fact that two arbitrary nodes can be con- 
nected through a relatively small number of steps, known as the small world effect. We calculate the 
probability that two arbitrary nodes in the network can successfully communicate with a fidelity 
above a given threshold. This amounts to working out the classical problem of percolation with 
limited path length. We find that this probability can be significant even for paths limited to few 
connections, and that the results for standard (unlimited) percolation are soon recovered if the path 
length exceeds by a finite amount the average path length, which in complex networks generally 
scales logarithmically with the size of the network. 

PACS numbers: 03.67.-a, 03.67.Bg, 89.75.Hc, 64.60.ah 



I. INTRODUCTION 

Networks permeate all informational structures. They 
underlie natural, social and artificial systems where dif- 
ferent parties interact, describing the flow of information 
between them. Differences in the characteristics of such 
interactions and how they evolve give growth to different 
types of structures: regular lattices, completely random 
networks or, spanning the range between these two, com- 
plex networks, which do not have a regular structure but 
neither are completely random. Quantum information is 
not an exception, and quantum networks [Tj where nodes 
communicate between them through quantum channels 
are essential to quantum information processing and dis- 
tributed applications. One of the key tasks in these net- 
works is the transmission of quantum information be- 
tween two distant nodes of the network. This task de- 
pends not only on the quality of the connections between 
nodes and on the amount of resources, but also on the 
underlying structure of the network. Therefore, under- 
standing how structural properties affect the functional- 
ity of the network will allow both the design of better net- 
work architectures and the modification of existing ones 
that make feasible communication at further distances, 
among a greater number of nodes or in the presence of 
higher levels of noise. 

Two distant nodes in a network may be connected by 
one path of entangled states (Figure [I]). In this case, 
long-distance entanglement between two nodes can be 
established with a probability that decays exponentially 



with the distance separating the nodes. This problem 
can be overcome by quantum repeaters, which create a 
distant entangled pair of high fidelity [2j. However, such 
technique require a number of qubits in each node that 
scales logarithmically with the distance [3] . More impor- 
tant, though, is that it only considers a one-dimensional 
connection between the two nodes. These two nodes, 
however, may be embedded in a more realistic, higher 
dimensional network. In this higher number of 

paths may exist which can help in the communication: 
with the existence of clusters of nodes connected by en- 
tangled states, two distant nodes will be able to estab- 
lish entanglement between them if they both belong to 
the same cluster [4]. Entanglement percolation, which 
makes use of such higher dimensional networks, was first 
proposed in the honeycomb lattice [5\ and later extended 
to other regular lattices [7], to schemes using multi- 
partite entanglement [8] and to noisy networks [9TH2] . In 
Ref. [13] , we studied entanglement distribution in a wide 
class of complex networks with pure state connections. 
Complex networks arise in many real scenarios, notably 
including the most important real world communication 
networks, and it is very plausible that they will become 
relevant in quantum communication architectures too. 
They offer very rich properties and phenomena. Interest- 
ing quantities can be computed requiring only statistical 
properties. This might seem a limitation but it can rep- 
resent an advantage: it makes mathematically tractable 
some problems that are hard or impossible to solve on 
lattices, and provides a minimal description in scenarios 
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FIG. 1. (Color online) A quantum network. Nodes are big 
(green) circles. Links between nodes are represented in grey, 
each holding a number of bipartite entangled states (black 
lines) shared between nodes. 



were complete knowledge of the system is not available 
or is hard to obtain. 

One of the primary features of networks is the pres- 
ence, or absence, of a cluster of nodes connected between 
them (a connected component, in graph theory language) 
whose size is of the order of the size of the network. Such 
cluster is called the giant connected component, and in 
the asymptotic limit of infinite size networks it is de- 
fined as the cluster spanning a finite fraction S of nodes 
of the network. This concept is very closely related to 
that of a percolating cluster. In bond percolation, for 
example, edges are occupied with some probability </>i, 
thus connecting their end vertices, and empty with prob- 
ability 1 — <f>x, disconnecting them. Then there exist a 
percolation threshold 4>l in this occupation probability: 
below the threshold, all components are of finite size, 
while above it there exist one giant connected compo- 
nent whose size is comparable to the network size. This 
threshold is the critical point of a phase transition, gen- 
erally of second order, and can be manifested by the di- 
vergence of the average component size, which acts like 
a susceptibility in a magnetic material. 

The percolation threshold and the size of the giant 
connected component, as well as many other properties, 
strongly depend on the basic structure of the network 
[T3HH] as well as on degree-degree correlations [IT] [TS] 
and clustering |19j . Therefore, a change in the structure 
of a network can affect its ability to communicate infor- 
mation. For example, the scale free topology of Internet 
makes it strong — resilient — against the failure of random 
nodes [20], but not against target attacks directed to its 
major hubs (nodes with highest number of neighbors) 
[2"Tl |2"2"] . This relation between the structure of the net- 



work and the communication over it can be also exploited 
to benefit the earlier appearance of the giant cluster and 
to find architectures that allow communication even in 
the presence of noise. 

In this paper we study the distribution of quantum in- 
formation over quantum complex networks. We first fo- 
cus on networks where nodes are connected by bipartite 
pure entangled states. We propose a transformation of 
the network that, using only local knowledge, can change 
the structure of the network and decrease its percolation 
threshold. We also calculate how the percolation thresh- 
old and the size of the giant component change after the 
transformation. Then, we turn to mixed state connec- 
tions between nodes and show that the small world be- 
havior of many complex networks allows quantum com- 
munication above some fidelity bound for finite, but very 
large, quantum complex networks. 



II. RANDOM GRAPHS 

A network is naturally represented by a graph G, which 
is an ordered pair of sets G = {V, E}. V is the set of ver- 
tices (or nodes, or points), and E the set of edges (or 
links, or lines), which are pairs of elements of V and rep- 
resent the connections between them. In this paper we 
consider always undirected graphs with neither multiple 
edges (i.e., either zero or one edge between every pair of 
vertices) nor self-loops. The degree of a vertex, k, is the 
number of edges emerging from it. A connected compo- 
nent, or cluster, is a subgraph where any two vertices are 
connected by at least one path of edges and to which no 
more vertices can be added without losing this property. 

Random graphs [22] are ensembles of graphs Q of the 
same size, with a probability P(G) assigned to every 
graph G in the ensemble. For any property 0{G) of a 
graph we can calculate its average over the ensemble: 

(0) G = ]T 0(G)P(G). (1) 
GeQ 

However, in most real scenarios, but also in many theo- 
retical models, only a single, large graph is studied. Such 
graphs are said to be self-averaging if the property we 
are studying is well characterized by its mean. This hap- 
pens when the graph is large enough to make fluctuations 
around the average vanish. For a more detailed discus- 
sion about self- averaging, see e.g. Ref. [53] for random 
graphs and [S3] for the World Wide Web network. We 
will consider this assumption in the following, and check 
its validity numerically in the examples we consider. 

One of the basic properties of a graph is the distri- 
bution of the probability pk that a vertex has degree 
k. There is also a related distribution that will come in 
handy later, and is that of the excess degree: the num- 
ber of edges, k, emerging from a vertex reached through 
another edge, and excluding it. This probability can be 
found easily by considering first the degree of a vertex 
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reached through an edge. Since vertices with higher de- 
gree are easier to reach, such probability is proportional 
to the degree of the reached vertex, kpk / {k} . The excess 
degree probability is therefore 



(k + l)p fcH 
(k) 



(2) 



Generating functions [26] are a mathematical tool that 
shows very useful when studying properties of graphs de- 
scribed by probability distributions [IS]. Among other 
useful properties, they allow for a straightforward con- 
volution of distributions. Let us introduce them using 
the degree distribution. The function g p {x) that gener- 
ates the distribution {pk} is the power series of x with 
coefficients equal to the probabilities in the distribution: 



9 P ( X ) = ^Pk3 



(3) 



fc>0 



Note that each probability p^ can be recovered from its 
generating function ([3| by taking the k-th derivative of 
g p {x) at x = 0, 



Pk 



d k g p (x) 



dx k 



(4) 



x=0 



Since the probability distribution is normalized, ^2 k Pk = 
1, so is its generating function, g p (l) = 1. It is also 
convergent for \x\ < 1, which is all what we will use here. 

The first moment of the distribution p^ , which corre- 
sponds to the average degree of the graph (k), is equal 
to the first derivative at x = 1: 



(5) 



k>l 



This allows to express the generating function for the 
excess degree distribution, g r (x), in terms of 



g r (x) 



^r k x k 



(G) 



Higher moments can be similarly found by taking more 
derivatives. In general, the n-th moment is 



(k n ) = k >k = 



k>0 



Ax 



g P (x) 



(7) 



Convolution of independent distributions can be ob- 
tained by multiplication of their respective generating 
functions. For example, the total number of edges emerg- 
ing from n independent vertices (the sum of their degrees) 
is generated by [g p (x)] n . 



III. NETWORK EXAMPLES 

In this paper we calculate properties such as the av- 
erage component size, the giant connected component 
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FIG. 2. (Color online) Example of the first five crowns of a 
Bethe lattice with coordination number k = 3. Left: Bethe 
lattice before g-swap. Right: Same network after applying 
3-swap. 



size and the percolation threshold. Analytical results are 
found for random networks with uncorrelated degree dis- 
tribution, and for the Watts-Strogatz small world model 
in the mixed state scenario. We also discuss several net- 
work topologies as concrete examples of our results: the 
simple Bethe lattice, two networks (Erdos-Renyi and 
scale free) belonging to the configuration model, the 
Watts-Strogatz small world model and two real world 
networks. All of them share a common property known 
as the "small world effect": the average path length, or 
intervertex distance, scales logarithmically with the size 
of the network, rather than as a positive power of the size, 
N 1 ^, as is the case in finite-dimensional networks. Here 
we present a short description of each of these network 
models. 

Let us start with the simplest network example. The 
Bethe lattice is not a random graph but has some simi- 
lar properties such as a local tree-like structure and the 
small world effect (as long as the degree of its vertices 
exceeds 2), while at the same time remains amenable to 
analytical study. A Bethe lattice with coordination num- 
ber k is defined as an infinite regular graph where every 
vertex has the same degree k and is topologically equiv- 
alent to all the others, as shown in Figure [2j Random 
regular graphs — graphs where all vertices have a fixed 
degree but edges are placed randomly — asymptotically 
approach Bethe lattices, making them a relevant model 
where analytical treatment is usually possible. 

Erdos-Renyi graphs [27H29] are maximally random 
graphs with the only constrain (k) = z. An Erdos- 
Renyi network with N vertices can be realized by ran- 
domly placing M — Nz/2 edges, or similarly by placing 
an edge between every pair of vertices with probability 
z/N (also known as the Gilbert model), which is asymp- 
totically equivalent [2>] . Figure [3] shows an example of 
a small Erdos-Renyi network. Their degree distribution 
is Poissonian, p k = e~ z z k /kl, with generating functions 
g p {x) = g r (x) = exp[z(x - I)]. 

Real world networks are not Poissonian but typically 
exhibit a power-law (scale free) degree distribution, p k ~ 
k~ T , characterized by a relatively important number of 
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FIG. 3. (Color online) Left: Example of an Erdos-Renyi 
network with N = 30 and z = 4/3. Right: Example of a 
Watts-Strogatz network with N = 30 and P = 0.2. 



nodes with a degree much greater than the average. Scale 
free networks with r < 3 have a percolation threshold at 
4>i = 0, while for networks with r > 3 a finite thresh- 
old appears. However, in heavy-tailed networks such like 
these, a cutoff in the degree naturally appears in scenar- 
ios where high degrees cannot exist due to, e.g., targeted 
attacks, physical constraints, saturation effects or finite 
size networks. For this reason we consider scale free net- 
works with an exponential cutoff, pk = Ck~ T e~ k / K (C is 
a normalizing constant) , while the pure scale free behav- 
ior can still be recovered by taking the limit k — > oo. The 
cutoff k strongly affects the network properties, and in 
particular networks with r < 3 have now a finite thresh- 
old. 

Random graphs with uncorrelated degree distribution 
such as the previous models exhibit a very low level of 
clustering (also known as transitivity): the likelihood 
that two neighbors of the same vertex are also neigh- 
bors between them. Aside from regular lattices, which 
have a high level of clustering, there are also many real 
world networks with this property. This is specially true 
for social networks, but also for communication and bi- 
ological networks. To study this behavior, many mod- 
els have been proposed, maybe the most studied being 
the Watts-Strogatz model |30j . This model is a random 
graph with ordered local structure and high level of clus- 
tering but still with surprisingly low average path length 
(see Figure [3j Here we study a slight modification, also 
considered in [14]. It is generated by placing N vertices 
in a ID ring. Then, N additional random edges, called 
"shortcuts" , are added with probability /3, thus giving an 
average of j3N shortcuts. 

We also study two real world networks. The first case 
represents a real world scale free network consisting of 
World Wide Web sites in the nd.edu domain [31]. In this 
case we introduce an artificial cutoff by neglecting nodes 
with degree k > 15, leaving a graph with 442 192 nodes 
and 170 352 edges. The second real world example is the 
OpenPGP Web of Trust, a social network representing 
the trust between OpenPGP users. Without going into 
much detail, OpenPGP is a standard encryption proto- 
col for securing email communications using public key 



cryptography. If Alice wants to send a secure message to 
Bob, she has to use Bob's public key to encrypt it. The 
authentication problem arises when Alice cannot verify 
if the key she is using is really owned by Bob. A solu- 
tion to this problem is the Web of Trust, in which every 
user signs a public key if she trusts it, thus generating 
a directed graph. To trust a key, usually a user has to 
meet with the key owner and check that he is really who 
he claims to be. This social model is thus relevant to 
quantum communication in the sense that at this point 
the two users could create a bipartite entangled state 
and then separate, each keeping one of the parts. By 
repeatedly doing so between different pairs of users, as 
in the Web of Trust, a quantum network would be cre- 
ated. Here we use the strongly connected component of 
the Web of Trust obtained from the Swiss keyserver [? ] 
as of May 25, 2010, containing 41 459 keys and 424 577 
signatures. We considered only bidirectional edges, cor- 
responding to users who mutually signed their keys. This 
leaves an undirected graph with 38 550 keys and 145 388 
two-way signatures. 



IV. PURE STATE NETWORKS 

We first focus on networks of pure, nonmaximally en- 
tangled states, as in [5]. In this case, edges have some 
probability <f> n of being converted into maximally entan- 
gled states depending on the amount of entanglement 
and the number n of bipartite states per edge. Singlets 
can then be used for perfect teleportation, i.e., they are 
equivalent to a single-use ideal quantum channel. This 
strategy can be directly mapped into a bond percolation 
problem, and is thus called classical entanglement perco- 
lation. There exist then a critical probability <j>\, called 
the percolation threshold, above which the giant com- 
ponent appears with fractional size S > 0. Above this 
threshold, any two nodes are able to share maximal en- 
tanglement if they both belong to the giant component. 
This happens with a probability S 2 which is indepen- 
dent of the distance, but that strongly depends on the 
network topology. Hence, for a given type of edges 4>„ 
long-distance entanglement will only be possible for net- 
works fulfilling (j> n > <fi*. Remarkably, due to the quan- 
tum nature of the connections it is possible to drastically 
change the network topology by local actions: a particu- 
lar measurement is done on qubits within the same node, 
establishing new connections between neighboring nodes. 
Thus, a quantum preprocessing of the network can be 
carried before edges are converted into singlets, so the 
new structure provides, e.g., a better percolation thresh- 
old. Moreover, in order to carry the particular prepro- 
cessing strategies we consider it is not necessary to know 
the precise structure of the network. Given only general 
statistical properties of the network, we propose strate- 
gies that act on each node depending only on locally ac- 
cessible information, such as the degree of the node. We 
calculate the new percolation threshold (f>t of the mod- 



ified network and the evolution of the giant connected 
component S. Above the threshold any two nodes will 
be able to establish an entangled state with probability 
S 2 , again independently of the distance between them. 
Thus, quantum preprocessing can benefit communication 
in two ways: by lowering the percolation threshold and 
by an increase in the giant component size. 



A. Network model 

We consider a quantum network in which neighboring 
nodes share n — 2 copies of a bipartite pure entangled 
state of two qubits, 



|^) = VAo|00) + VAi|ll>, 



(8) 



where ^/Ao > \/Ai > are its Schmidt coefficients. A 
partially entangled state can be converted into a maxi- 
mally entangled state (singlet for short) with singlet con- 
version probability (SCP) that only depends on its largest 
Schmidt coefficient 132 . For the state lib), the SCP is 



^=mm[l,2(l-Ao)]. 

We will consider edges that are of the form 
thus can be converted to singlets with SCP 

02 = min[l,2(l-Ag)]. 



(9) 



and 



(10) 



With this probability two neighbors can establish a per- 
fect channel between them. As we discussed above, for 
two distant nodes this probability depends on the struc- 
ture of the network that connects them. 



B. Modifying the network: g-swap 

In |13j we introduced a network transformation, the 
q-swap, that requires only local information of the net- 
work: the degree of a target node and the status of its 
neighbors. The g-swap is built upon a basic transforma- 
tion: entanglement swapping or swap 33] . In a subgraph 
with three nodes, the party at the target node c performs 
a Bell measurement on two qubits, each of them belong- 
ing to states shared with different nodes a and b (see 
Figure [I| . After this operation, the central qubits be- 
come disentangled from a and b, but in return a mixed 
entangled state with the same SCP as \ip) is created be- 
tween a and b [S]. Note that this operation can not be 
repeated with a fourth node because the newborn state 
shared between a and b is not of the form of The 
q-swap performs swap transformations between succes- 
sive pairs of neighbors of a central target node of degree 
q, thus changing and initial g-star with edges {ip)® to 
a g-cycle with newborn edges, while the central target 
node becomes disconnected from the network (see Fig- 
ure El) . For a given network topology, we will see that 




FIG. 4. (Color online) Entanglement swapping (a) and q- 
swap (b). Black dots represent qubits, big circles nodes, solid 
lines states \ip) and dashed lines the resulting mixed state 
after entanglement swapping, with same SCP as \ip). 



performing g-swaps on nodes with certain degrees im- 
proves the threshold. It is worth noting, however, that 
in some instances the application of particular q-swaps 
may be counterproductive. 



C. Percolation threshold and giant component 

The main two figures of merit that we will use to com- 
pare the two strategies are the percolation thresholds ((f)* 
and <f>l for the classic and g-swap strategies respectively) 
and the size of the giant components (S and S). Since q- 
swaps disconnect vertices, which can be chosen not to be 
the two corresponding to the parties that want to com- 
municate, the probability to connect two remote nodes 
is in fact S* 2 , where S = SSi/Si and Si is the value of 
S at <pi = 1. The percolation threshold tells us which is 
the minimum amount of entanglement needed for long- 
distance communication, while the square of the giant 
connected component size is the probability that any two 
nodes can communicate. 

To compute these two values we will use the generat- 
ing function formalism described in Section [TTJ. The key 
probability distributions are those of finding a connected 
component of finite size s, either when a random vertex 
is chosen, P s , or when a random edge is followed to one 
of its ends, R s . A random edge is empty with proba- 
bility i?o = 1 — 02! giving a cluster of size 0. When it 
is occupied, then a node of degree k + 1 is reached with 
probability r^, giving access to k clusters. Therefore, 

oo 

k=0 si)S2,-' ,sjj 

(11) 

We have assumed that components are treelike, i.e., that 
they do not have finite loops. This is indeed true for finite 
components, since an edge exiting such a component will 



reconnect back to itself with probability proportional to 
s/N —> 0. The function generating R s is then h R (x) — 
J2 s >o RsX s ') which gives the recurrence relation 



h R {x) = 1 - 02 + 4>2xg r [hR{x)}. 



(12) 



It is crucial to notice here that by restricting to finite 
s we have explicitly excluded the infinite giant compo- 
nent from hp(x). Thanks to this the previous treelike 
assumption holds. 

We can proceed similarly with P s . A random vertex 
has degree k with probability pk, giving acces to k clus- 
ters. The probability that this random vertex is in a 
component of size s is then 



P, 



Pk 



E 



R Sl R s 



(13) 



Now the generating function for P s is related to hp_(x), 



h P (x) = xg p [h R (x)}. 



(14) 



Knowledge of hp{x) and h R (x) allows for the deriva- 
tion of 4>l and S. The probability u that an edge con- 
nects to a finite component is the smallest real solution 



of u = h R (l) in Eq. (12 1, which is in general a transcen- 
dental function. In fact, the percolation threshold is the 
value of 4>i at which a solution u < 1 appears. More- 
over, u 2 is the probability that the edge connects to a 
finite components through both ends, so with probabil- 
ity 1 — u 2 a random edge belongs to the giant component 
(which is below the threshold). 

Similarly, the main quantity of interest, the giant com- 
ponent size, can be computed as the missing component 
in the whole network. The sum of all P s , hp(l), gives 
the probability that a random vertex belongs to a finite 
component, while the probability S that it is in the giant 
component is S = 1 — J2 s >i Ps = 1 — hp(l). Again, this 
equation is usually transcendental and has to be solved 
numerically. For instance, the Erdos-Renyi model, with 
g p {x) — g r (x) = e ztyX ~ x \ has a giant component fraction 
S = 1 — e~ z ^ 2S [15]. In this case the solution can be 
expressed in terms of the Lambert W function, 



S = 1 + —W(-zfae 

Z(f>2 



-Z<t>2 



(15) 



and the phase transition to S > occurs at the well- 
known point 4>\ = 1/z. On the other hand, first moments 
can usually be computed even when a closed expression 
for hp(x) and h R {x) is not known. As a relevant example, 
the average component size (s) = h' P (l) is an important 
property of the network that provides an alternative way 
of finding the probability threshold: it is at this point 
that (s) diverges. From the derivatives of Eqs. (12 1 and 



( 14 1 it is immediate to find that this divergence can be 



traced back to that of 




(b) 



•X 




FIG. 5. (Color online) Example of the branching process be- 
fore and after a 3-swap, starting at the leftmost node, (a) 
Before any operation the branching process arrives at a node 
of degree 3, leading to 2 components (in dark grey), (b) After 
the 3-swap, the branching process is already in a 3-cycle, each 
of its nodes belonging to one of the 2 components. 



This brings the general result for the critical SCP 4>* 2 — 

i/^(i) m- 

We now want to understand how the q-swap trans- 
formation changes the percolation properties of the net- 
work. Every particular q-swap can be implemented (or 
not) with probability U q (or 1 — H q ) on nodes of degree 
q. Giving the values for each H q specifies the quantum 
strategy, q-swaps introduce cycles, so components are 
no longer treelike and generating functions can not be 
directly used. Note however that, since newborn edges 
cannot be reused, those cycles do not overlap between 
each other and can thus be treated as blocks of a tree- 
like component by considering two steps in the branch- 
ing process. We first compute the generating function for 
the probability R s after q-swaps are done, h R {x). Now, 
instead of arriving to a vertex of degree q connecting to 
other q— 1 components, after a q-swap operation has been 
done we arrive at a cycle of q nodes (including the one 
we are coming from) connected via edges occupied with 
probability cj)i (see Figure [5]) . When edges are converted 
into singlets, the accessible nodes of this new q-cycle form 
a string of length / with probability 







(/ + i)#(i-&) s 



for I = q, 
for I = q 
for I < q 



1, 
2. 



For / < q — 2, I new components emerge, with total size 
(including all the vertices in the cycle, except the starting 
one) probability generated by [xg r (hp(x))] 1 . For I = q— 1 
and I = q, q — 1 components emerge, again with total 
size probability generated by [xg r (hp(x))) q ~ 1 . The total 
size of such cycle and its emerging components is then 
generated by 

q-2 



C q (x) = £(Z + 1)^(1 - 0i) 2 [xg r (h R (x)) 

1=0 

+ [qfi-Hl - <h) + <f>i] [x9r(h R (x))] 9_1 (17) 



(16) 



Therefore, the new hp(x) is of the same form of Eq. ( 12 ) 



7 




FIG. 6. (Color online) Two clusters of 5 nodes with degree 2. 
Nodes are big circles: empty if their degree is different from 
2, dark grey (red) if they are operated on, light grey (green) 
if they are not. Top: operations are done at nodes 1, 3 and 



5, leading to v^ 1 ^. Bottom: 
and 4, leading to Vi^ 1 ^ ■ 



operations are done at nodes 2 



plus a term h R ^ q (x) for each g-swap: 

h R (x) = 1 - <f> 2 + <t>2xg r (h R (x)) + ^ RqhR,q{x) (18) 



q>2 



h R , q (x) = T q - 



L) - <f> 2 x (h R (x)y 1 + C q (x) 



At this stage we can already calculate </>* as the smallest 
value of <pi for which there exists a positive solution u = 
h R (l) < 1 to ( 18 1 at x = 1. It is easy to convince oneself 



that each separate contribution h R ^ q {l) either increases 
or lowers the percolation threshold and therefore for the 
optimal strategy each H q is either or 1. 

For the new hp{x) we need to consider that not all 
nodes of degree q are suitable targets of g-swaps, since 
they cannot be performed on adjacent nodes. Therefore, 
given a node of degree q there is a probability rj q that a q- 
swap can be performed on it. If the g-swap is performed 
on a node, then it changes its degree from q to zero and 
hence 



h P {x) = xg p [h R (x)] + x^U q i] q p q {l 

q>2 



[h R (x)]*}. (19) 



By using the solution u = h R (l) here, we can obtain the 
size of the giant connected component, S = 1 — hp(l). 
This gives the probability S 2 that two distant nodes are 
connected by a path of singlets. 

The probability r\ q depends on which degrees are tar- 
gets of g-swaps and on how the network is traversed to 
operate on the nodes. To compute its value we need to 
consider maximal clusters consisting of nodes where all 
vertices are of any target degree q — the border of such 
clusters is necessarily made of nodes of degree different 
from g, and hence operations can be done independently 
on every cluster. As an example, let us discuss the sim- 
plest case of only performing 2-swaps. Starting from 
a random vertex of degree 2, we find a cluster of ver- 
tices of same degree 2 whose size is s with probability 
s(l — ri) 2 rj _1 . By acting on a node, and then on every 
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FIG. 7. (Color online) Probability r/2 of performing a 2-swap, 
given a verte x of d egre e 2. Upper and lower lines correspond 
to Eqs. (20 1 and (21 1 respectively, squares to Erdos-Renyi 



network simulations, circles to scale free network simulations. 
All the simulations were performed with N = 10 6 and Bread 
First Search traversal of the graph. 



second node, there are two possible values for the number 
of operations done in each cluster, \s/2] and [s/2\ (Fig- 
ure [6]) , which coincide for s even. This gives a maximum 
and minimum value for 772, 



(max) 



-1 \s/2] 



1 



(l-ri) 3 + 5>(l-ri) 2 r; 



1 + 7-1 

-1 [s/2\ 



(20) 



s>2 
2 



1 - (1 - r Qrj 

1 + 7-1 



(21) 



Note that for clusters of size s = 1, an operation is always 
done. When operations are performed starting from a 
random vertex in each cluster of vertices with degree 2, 
one needs to take into account the number of vertices s 
and t at odd and even (including zero) distance from the 
first vertex: operations will be performed on a fraction 
t/(t + s) of the cluster. The probability £(s, t) of starting 
in a vertex of degree 2 such that it has s neighbors of 
degree 2 at odd distance and t at even distance is 



£(M) = 



1 



2 

■s-t 



(22) 



if \s — 1\ < 1 and otherwise. For general q, this proba- 
bility can be found by generating functions similar to the 
ones described before, see Appendix |A"| for more details. 
Given the probability £(s,i), then the value for 772 when 
operations are started at each cluster of target vertices is 



(rand) 



^ t- 

s,t 



7"l + (1 — ri) 2 atanh(ri) 
2^ ' 



(23) 

Figure jlshow ?7 2 max ' 1 , r^ 1 ^ and 772™"^ together with 
numerical simulations performing 2-swaps by traversing 




FIG. 8. (Color online) Erdos-Renyi network, (a) Gain 7 as a 
function of the mean degree z after 2-swap (dotted line), 2,3- 
swap (dashed line), 2,3,4-swap (dot-dashed line) and optimal 
g-swaps (solid line), (b) Normalized size S of the GCC as a 
function of <j>i for z — 2.5, before (squares) and after (circles) 
2,3-swap, N = 10 6 . 




FIG. 9. (Color online) Scale free network, r = 1. (a) Gain 7 
as a function of the cutoff k after 2-swap (dotted line), 2,3- 
swap (dashed line), 2,3,4-swap (dot-dashed line) and optimal 
q-swaps (solid line), (b) Normalized size S of the GCC as a 
function of </>i for n — 4, before (squares) and after (circles) 
2,3-swap, N = 10 6 . 



the graph with a Breadth First Search, as described in 
the following section. The numerical values for 7? 2 are 
close to the maximum value because it is much more 
likely that the traversal of graph started outside most of 
the degree 2 clusters (e.g., arriving through one of the 
white nodes in Figure [6]) , thus performing the maximum 
number of operations in them. 



D. Network examples and simulations 

Here we present some examples of entanglement per- 
colation in the networks described in Section |III| and 
we provide analytic solutions for paradigmatic cases. To 
check these results and extend them to correlated and 
real world networks we have performed computer simu- 
lations with various networks models. Graphs with un- 
corrected degree distribution are relatively easy to gen- 
erate [15]. First, a set of N numbers {k t } randomly cho- 
sen to follow the desired degree distribution is generated, 
so each vertex i has ki stubs or "half edges" associated 
with it. If the sum ki is odd, a new set is gener- 
ated until an even sum is obtained, so all stubs can be 
joined. Then pairs of stubs are selected randomly and 
joined to form edges until there are no stubs left. In our 
simulations we did not allow self-loops or multiple edges. 
The quantum preprocessing is done by traversing all the 
graph with a Breadth First Search (BFS), which starts 
at a random root vertex and explores all the neighboring 
nodes at distance 1, 2. . . , in order, until all the vertices 
in the component have been visited. After that, another 
BFS is done starting from a random unexplored vertex 
in another component, until all components have been 
examined. At each discovered vertex, the local structure 
is changed from a g-star to a g-cycle if the vertex degree 
is one of the target degrees and if non of the edges in the 
star have already been used. 

For the Bethe lattice with coordination number q, 
g p (x) — x q and g r {x) — x q ~~ x . In this network the phase 
transition occurs at (f>2 = {Q ~ 1) and 

(q - l)- 1 = (1 - 0i)- 1 {2^i + 4>l[MQ - 1) - (« + 1)]} 



after g-swap is applied. Therefore, g-swap gives always a 
better threshold except for the special case q = 2 of an 
infinite ID chain, where the probability decays exponen- 
tially with the distance. 

In the Erdos-Renyi network, before any transforma- 
tion the threshold is given by = 1/z. After, e.g., the 
2-swap and 3-swap operations, the thresholds are, respec- 
tively, 



4>2 



-0 2 + z(20i -<j>t)} 



(24) 



and 



- -0 2 + ze- 2 h^ 2 + z(l + <h. -<%)]■ (25) 

Figure [8] shows the evolution of the giant connected com- 
ponent before and after the transformations, with perfect 
agreement between analytical and numerical results. Fig- 
ure ^ also shows the gain 7 = (0* — <f>i)/<fi* in the perco- 
lation threshold, which in some situations is higher than 
20%. The performance of different g-swaps depends on 
the mean degree z, usually improving the threshold those 
operations which act on nodes whose degree is around z. 
Figure|9]show similar results for the giant connected com- 
ponent evolution and the gain in scale free networks with 
r = 1. In this case the gain can be of around 25%. 

For the Watts-Strogatz model, which is correlated, and 
the World Wide Web network, the above approach is 
not valid because the tree-like assumption does not hold. 
However, numerical simulations show that g-swaps can 
also provide an improvement in the percolation thresh- 
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Figure 



10 



shows the threshold probability 
for the Watts-Strogatz before and after 2-swap and the 
size of the giant connected component. Figure 11 shows 
the size of the giant connected component for the World 
Wide Web. 

Note that, in general, it may be counterproductive to 
perform q-swaps. In the above figures we see that for 
some values of <f>± the giant connected component frac- 
tion S without preprocesing is larger than S. This of- 
ten happens around 4>\ = 2 — ^/2. This is precisely the 
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FIG. 10. (Color online) Watts-Strogatz network, (a) Per- 
colation threshold cj)* as a function of shortcut probability /3 
before (squares) and after (circles) 2-swap. The solid line is 
the analytic result from Ref. |14] , (b) Normalized size S of 
the GCC as a function of 0i for /3 = 0.2, before (squares) and 
after (circles) 2-swap, N = 10 s . 




FIG. 11. (Color online) World Wide Web network [31], with a 
cutoff at k — 15, before (squares) and after (circles) 2,3-swap. 



point where the edges in the unmodified network can be 
directly converted into singlets with $2 = 1, i.e. all con- 
nections become ideal channels and S attains its maximal 
value S = S\. Obviously at this stage any preprocessing 
cannot further increase the size of the connected compo- 
nent, and it will most likely decrease it. 



V. MIXED STATE NETWORKS 

A network with nodes connected by pure states is 
an abstraction that gives insight into the possibilities 
of long-distance entanglement in complex networks, en- 
abling perfect teleportation between distant parties when 
at least a path of maximally entangled states is created. 
In general, however, states connecting two neighbors are 
noisy and need to be described by mixed states. In this 
situation, the optimal fidelity of teleportation / is di- 
rectly related to the maximal singlet fraction F [35] by 



/ 



Fd+1 
d+1 



(26) 



where d is the dimension of each part of the bipartite 
state and F is defined as the maximal overlap of a state 
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FIG. 12. (Color online) For a limited path length, cluster 
growth depends on the network topology. Here, clusters of 
limited path length I = 1,2,3,4 in the OpenPGP Web of 
Trust, with sizes 11, 37, 115, 286. 



p with a maximally entangled state |\&) = ^ J^i 



F(p) = max<*|p|*} 
I*) 



(27) 



Entanglement percolation in the mixed state scenario has 
already been addressed in regular lattices with connec- 
tions consisting singlet that have suffered an amplitude 
damping [9j[l0]. There, an hybrid swapping strategy is 
proposed for this type of connections, which could also 
be used to build a mixed state g-swap to act on com- 
plex networks with at least four states per edge. In an- 
other approach, Perseguers gives a fidelity threshold for 
the links above which long-distance quantum communi- 
cation in the presence of noise is possible for an infinite 
cubic lattice [Hj. 

Communication in noisy networks can be considered 
from another perspective. The noise in the connections 
fixes a limit I in the maximum number of nodes through 
which the information can be repeated before it becomes 
too corrupted In this limited path length scenario, 
the total number of vertices that a given node can com- 
municate to also depends strongly on the structure of the 
communication network (Figure 12). This is related to 
the average path length Z av : the length of the shortest 
path averaged over all possible pairs of nodes. All nodes 
within this distance constitute a significant fraction of 
the network. Therefore, for a path length limit I above 
the average / av , communication will be possible among an 
important number of nodes. Since the limiting I is finite, 
the giant connected component appears only in models 
where l av is also finite. In general, this only happens if the 
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network size is finite too. The question then is whether 
a small I will suffice to cover a significant fraction of the 
network. In finite <i-dimensional networks, the average 
path length scales as l av ~ N 1 ^. However, the average 
path length of many complex networks scales logarithmi- 
cally with the size of the network. This property is known 
as the small world effect, and appears also in many real 
world communication networks such as Internet. In this 
case, to access a significant fraction of nodes, only a small 
number of edges need to be traversed. Small world mod- 
els are therefore the first candidates, where loses by noise 
can be balanced by a short path length. 

The problem of limited path percolation was also ad- 
dressed in a different approach by Lopez et al. [3B]. In 
their model, they calculate the percolation phase tran- 
sition under the assumption that communication is only 
effective if the new minimum path length between two 
nodes do not exceed a multiple of the original path length 
between them. Thus, in their study the limitation in the 
path length comes from the topology of the network and 
not from the nature of channels connecting nodes, which 
fixes a constant limit of nodes through which the infor- 
mation can be repeated. 

Here, we are interested in the number of nodes that 
can exchange quantum information with a given node 
for some fixed minimum fidelity, or similarly with what 
probability two random nodes can reliably communicate 
between them. We will consider a similar scenario as 
in the previous sections, but replacing pure-state con- 
nections with generic entangled mixed states. Here, no 
quantum preprocessing will be possible. However, we 
will find that the complex network structure (in par- 
ticular the small world effect) allows to interconnect a 
large number of nodes using the standard entanglement 
percolation strategy. We start by doing some numerical 
simulations and then derive the generating functions for 
limited path percolation and compute the limited average 
size in non-correlated networks and the Watts-Strogatz 
model. 



A. Network simulations 

We begin by simulating different models of networks. 
For simplicity we consider that edges hold a single copy 
of a two qubits state pp with maximum singlet fidelity 
F > 1/2 so that the classical limit of / = 2/3 in the tele- 
portation fidelity can be exceeded. Note that, as long 
as pf is entangled, this limit can be achieved even if 
F < 1/2 by locally increasing the singlet fidelity through 
trace-preserving local operations and classical communi- 
cation (LOCC) [37 . By applying random bilateral rota- 
tions, pf can be brought into a Werner state, 
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FIG. 13. (Color online) Normalized Z-limited average compo- 
nent size {si)/N as a function of I — i av for the Erdos-Renyi 
network with k = 2 and network sizes N = 10 3 , 10 4 , 10 5 , 10 6 
(squares, circles, diamonds and triangles). Superposed filled 
markers (in orange) is the shape of the path length distri- 
bution normalized with the total number of possible vertex 
pairs, n(l) = 2n(l)/N(N - 1). Solid black line is Eq. pi) , 
horizontal dashed line is the square of the GCC at <f>i — 1, 



see Eq. ( 15 1 



which has the same singlet fidelity F. This state can also 
be written as 

p a =a|¥-)(¥-| + (l-a)i (29) 

with a = (4F— 1) /3. It can be interpreted as the result of 
transmitting a pure singlet through a depolarizing chan- 
nel. Hence, when a state is teleported through I of such 
edges [H O [35] , hs fidelity is 1 with probability a 1 and 
1/2 otherwise, so its final fidelity is fi = (1 + cr)/2. This 
fidelity decreases exponentially with the distance I and 
makes such communication scheme useless in networks 
such as linear chains or regular lattices, where the typi- 
cal distance between two nodes scales as the size of the 
network. However, as we discussed, the typical distance 
in many complex networks scales only logarithmically. 
The maximum distance I that information can travel is 
fixed by the minimum fidelity / m i n required at the end 
point and by the purity a of the channels, 



/ = 



ln(2/ min - 1) 
In a 



(30) 



This means that, even if there exists a path between a 
sender and a receiver in a network, it will only be useful 
if the length of this path is below a certain threshold. 
We performed extensive simulations of networks where 



neighboring nodes share a state ( 29 1 , and considered the 



(28) 



classical limit as the minimum required fidelity, / m i n = 
2/3. For small networks (N < 10 4 ) we performed the cal- 
culations over several network realizations and then av- 
eraged the results. For bigger networks, a single network 
realization is usually enough due to the self-averaging. 
The Z-limited average cluster size (si) is a specially rel- 
evant parameter, which amounts to the probability that 
two nodes can communicate with fidelity / > / m i n . 
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FIG. 14. (Color online) Normalized i-limited average com- 
ponent size (si)/N as a function of (I — l^N^ 1 ^ 2 for the 
Honeycomb 2-dimensional network with network sizes N = 
1014,5046,10086 (squares, circles and diamonds). Super- 
posed filled markers (in orange) is the normalized histogram 
of the path length distribution. 



N 




FIG. 15. Normalized /-limited average component size {si)/N 
as a function of singlet fidelity F in the biggest compo- 
nent OpenPGP Web of Trust (squares) and a honeycomb 
2-dimensional lattice (crosses), N ~ 3 ■ 10 4 . 



We have thus calculated (si) for different network mod- 
els and sizes. In Figure [13] we plot the normalized size 
(si) /N as a function of I — l av for the Erdos-Renyi model, 
with average path length l av ~ In Nj In z [39] . For differ- 
ent network sizes the curves collapse, supporting a linear 
TV-dependence (si) ~ N for fixed /. Similar results have 
been recently found for the average number of nodes at 
exact distance I from a random central node [40] . Re- 
garding the dependence in I, our results show that the 
average size grows exponentially with I for I<I JV , but 
deviate from this behavior when I is close to the aver- 
age path length, saturating to the maximum component 
size shortly after l av . This deviation is due to the deple- 
tion of nodes at distance I > l av . In the same Figure 13 



we plot the path length distribution, i.e., the number of 
pairs n(l) separated by a distance I, normalized by the 
total number of pairs N(N — l)/2. Again, both curves 
N = 10 3 and N = 10 4 collapse, thus supporting a de- 
pendence n(l) ~ TV 2 . We also found similar results for 
the scale free and the Watts-Strogatz models. It is in- 
teresting to note that, while 2* av grows with the size of 
the network, the width of the path length distribution 
remains constant. Thus, for large networks a small in- 
crease in a near l av leads to an abrupt change in (si/N). 
This is in stark contrast to regular lattices, where both 
the mean and the width scale as N x / d . For instance, in 
Figure[l4]we plot (si) /N and the path length distribution 
of the Honeycomb 2-dimensional lattice as a function of 
(I - l^N- 1 / 2 . The collapse of the curves confirms the 
N 1 ^ length-scale dependence. 

As an example of a real work network, we considered 
the OpenPGP Web of Trust. Figure [15] shows the prob- 
ability that two arbitrary nodes can communicate with 
fidelity / > 2/3 as a function of the singlet fraction F. 
Again, the comparison with a Honeycomb lattice of the 
same size shows that the small world property of the com- 
plex networks allows for faithful communication between 



most of the nodes in the network for reasonable values 
of the noise, while in regular lattices this is only possible 
for nearly pure states. 



B. Average component size in limited path 
percolation 



We now proceed to derive the generating functions for 
the limited path percolation problem. In this case, we are 
interested in the distribution of sizes s of the components 
that can be reached by only I steps through edges that are 
always occupied. As in the non limited case, there are two 
different distributions P^p and Rs for the cases where 
a random vertex or a random edge are selected. The two 
corresponding generating functions, h^) and , read as 




(2-1) 



(,;) 



and 



h<g(x) 




for I = 0, 
for I > 1, 



for I = 0, 
for I > 1. 



(31) 



(32) 



Note that all edges are occupied with probability one. 
The generalization to a different occupancy probability 
is straightforward, but not needed here. 

As before, we are now ready to obtain the Z-limitcd 
average size, 



{si) = 



dx 



= l + g' p (l)h't 1 \l) 



(33) 



By solving the recurrence equation given by h!k (1) with 
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l-L 



4>i = 1, and the size of the giant component Si tends 
to the non-limited size S±. Therefore, for networks with 
the small world property, i.e. l av ~ log N, one can inter- 
connect with a threshold fidelity (say, the classical bench- 
mark / = 2/3) any arbitrary pair of nodes in the network 
provided that the singlet fraction of the edges scales as 
F = 1 — 0(1/ log N) with the size of the network, which 
is clearly less stringent than the analogous constrain for 
d-dimensional networks F = 1 — 0(N" 1 ^ d ). 

We also consider the Watts-Strogatz model presented 
in Section |III[ which has a base circular lattice of size 
N with f3N randomly added shortcuts. In this case the 
derivation of the probability that a random vertex be- 
longs to an Z-limited cluster of size s, Pi , and its gener- 
ating function hp(x) uses the formalism of "local clus- 
ters" introduced in [14] . This "local clusters" are clusters 
in the base lattice (without considering the shortcuts). 
For a given I, the "local cluster" is always of size 21 + 1. 
Then, a shortcut at distance A — 1 from the starting ver- 
tex leads to a (global) cluster of size s' with probability 
P^!~ X \ A random shortcut emerges from the starting 
vertex with probability 1/N, from a vertex at distance A 
with probability 2/N, and lies outside the local cluster 
with probability (N — 21 + 1)/N. Hence, that shortcut 
will lead to a cluster of size s with a probability given by 
the generating function: 



FIG. 16. (Color online) Normalized i-limited average compo- 
nent size {si)/N as a function of I — / av for the Erdos-Renyi 
network (a) and scale free network (b). Points are simulation 
results for network sizes N = 10 3 , 10 4 , 10 5 , 10 6 (squares, cir- 



cles, diamonds, triangles) , solid lines correspond to Eq. ( 34 1 
horizontal dashed lines are the values of Si . 



the boundary condition /in. (1) = 1 we find 



W = 



1 



+ a' m 1 -^ 1 )) 

+ l-g' r (l) 



for I = 0, 
for I > 1. 



(34) 



This equals to the probability that any two nodes will 
be able to communicate with fidelity above / m ; n . Fig- 
ure 16 shows this result for the Erdos-Renyi and the 



scale free model, with very good agreement between the- 
oretical and numerical results below Z av . 

As we discussed above, this exponential growth of (si) 
is valid for I well below Z av . Our numerical simulations 
show that the validity of this approximation can be ex- 
tended to values near i av . Figure [T3| shows that the path 
length distribution is very peaked around Z av , and its 
width is independent of N. This implies, on one hand 
that our analytical approach holds true for values of I that 
fall out of this finite width (approaching from below) — 
see Fig. [16] On the other hand, the finite width implies 
that if £ is a few steps beyond ^ av then most of the nodes 
in the components will be reached before the limit dis- 
tance is attained. In this situation, Equations (31) and 




,('-!), 



X) 



2^h p 

A=2 



(X) 



There are 2/3N shortcut end-points that can similarly 
contribute to the total size of the cluster. Recalling that 
the generating function of the sum of sizes is the product 
of the generating function of each size, we find 



2/3N 



(35) 



where x^ l+1 ' is the generating function corresponding to 
the starting "local" cluster. In the limit of large N this 
can be simplified to 



l +2 ; -W[2l-l-h% X) (x)-2^2 l x=2 



(36) 

Again, we can obtain the limited average size by taking 
the first derivative at x = 1. For I = 0, (s ) = 1. For 
/ > 1, this results in the recurrence equation 



1-2 

2V(sa> 



1 + 2^ + 2/3 I (*,_!) 

+2 + 2/3 ««,_!> + <s,_ 2 )), 



A=0 



(37) 



which can be exactly solved. Figure 17 shows this result. 
We want to stress the fact that from these generating 
functions, (31 1 and (36), one can also calculate the prob- 



ability Ps up to any s by solving s + 1 iterations of them 



( 32 ) approach the non-limited case of ( 14 ) and ( 12 ) with and using Eq. @. 
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FIG. 17. (Color online) Normalized i-limited average com- 
ponent size (si)/N as a function of / — i av for the Watts- 
Strogatz network. Points are simulation results for network 
sizes N = 10 3 , 10 4 , 10 s , 10 6 (squares, circles, diamonds, trian- 
gles), lines correspond to the solution of Eq. (37 1. 



VI. CONCLUSIONS 

We have demonstrated that quantum complex net- 
works offer a powerful framework for entanglement dis- 
tribution in large systems. Regardless of their intricate 
structure, complex networks can be studied by their sta- 
tistical properties, which allows to analytically compute 
some interesting properties and to deal with them with- 
out knowing their exact structure. Here we have consid- 
ered entanglement percolation in networks where connec- 
tions are built on pure, non-maximal bipartite entangled 
states, and have studied a local quantum preprocessing 
of the network that can significantly decrease the perco- 
lation threshold and therefore allow quantum communi- 
cation for a lower level of entanglement. The quantum 
preprocessing we have proposed is local in two senses. 
First, quantum operations are done always on qubits that 
belong to the same node. Second, the decision whether 
or not to perform such operation depends on the local 
structure of the network (the degree of the target node 
and the status of its neighbors) and on information about 
general statistical properties of the network. We have 
calculated the percolation threshold, which marks the 
minimum level of entanglement needed to entangle two 
distant nodes with finite probability. We have also com- 
puted this probability, which amounts to the square of 
the giant connected component in the network. These re- 
sults are analytical for networks with uncorrelated degree 
distribution, and can be compared to previous results 
in classical networks, which shows that the preprocess- 
ing can substantially improve communication over such 
networks by manipulating its local structure. We have 
also studied numerically the Watts-Strogatz small world 
model and a real world network, and have found a similar 
behavior. 

In this approach, the links between nodes are pure 



quantum states. A more realistic scenario, however, 
needs to consider noise in the connections. Here we thus 
have considered the situation in which such connections 
are made of noisy mixed states. 

We have shown that in complex networks a direct im- 
plementation of the entanglement percolation strategy, 
without quantum preprocessing, allows for faithful quan- 
tum communication (above a fixed fidelity threshold) be- 
tween a large number of nodes. The noise severely limits 
the number of steps or connections through which infor- 
mation is transmitted. However, in complex networks, 
one can reach a sizable amount of nodes with a moder- 
ately low number of steps. If the fidelity threshold allows 
for a path length slightly higher than the average path 
length, all nodes in the giant component become faith- 
fully connected. The path length distribution is peaked 
at low values (scaling as log N in complex networks ver- 
sus TV 1 /'* in d-dimcnsional lattices), and has finite width 
(constant in N versus N 1 ^). This implies that in com- 
plex networks a finite fraction of faithfully connected 
nodes appears for much smaller limiting path lengths and 
reaches the giant component size abruptly. Hence, here 
the advantage of complex networks is twofold: the av- 
erage path length which marks the transition scales log- 
arithmically with the network size, and the additional 
steps needed to reach the non-limited scenario is finite. 

We have shown that new phenomena appear if net- 
works and the operations one can performed on them 
are governed by the laws of quantum mechanics. This 
has been known for regular lattices, but the rich proper- 
ties of complex networks still remain widely unexplored 
in the quantum setting. Our results in percolation, to- 
gether with new behavior found in the emergence of sub- 
graphs in quantum random networks [41) . are examples 
of these phenomena. 

Our results also contribute to the field of classical com- 
plex network. We have given analytical results for the 
gain in the percolation thresholds and the size of the gi- 
ant component for uncorrelated complex networks that 
undergo a set of local inversions (transformation that pro- 
duce the complement of the induced subgraph of the tar- 
get node) . The problem at hand of studying how critical 
properties of a network can be drastically modified by 
a given set of network transformations might be of gen- 
eral interest to other disciplines in the field. Finally, we 
have addressed the problem of limited-path percolation 
in uncorrelated and small world complex networks. 
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FIG. 18. (Color online) Left: Example of a connected com- 
ponent with three clusters (in dark grey) of degree 2 and 3. 
Right: Branching process in r/^ nd \ 3-swaps are made on 
dark grey (red) nodes, which are the t nodes at even distance 
from the top one. 



Appendix A: Calculation of r\ q 

As we said, the probability r\ q depends on the target 
degrees \q{\ and on how the network is traversed. By 
^(rand) we denote the probability rj q when a q-swap is 
first done in a random vertex with target degree, and 
then the cluster of vertices with degree belonging to {qi} 
is traversed by a Breadth First Search, performing q- 
swaps whenever possible (i.e., at every second step). Af- 
ter that, another vertex with target degree which has not 
yet been explored is selected, and its cluster traversed, 
until all target vertices have been checked. Such clusters 
consist of vertices of degree k £ {qi} that are connected 
by at least one path whose vertices have also a degree 
in {qi} and to which no more vertices of degree k can 
be added. Figure 18 shows an example of three of such 
clusters when the target degrees are 2 and 3. A ran- 
dom vertex of degree k € {qi} belongs to a cluster with t 
vertices at even distance (including itself) and s at odd 
distance with probability £(s,i). In this cluster of size 



s + t, t q-swaps are made. The probability 7 y^ rand ) j s then 



(rand) 



E 

t.s 



S + t 



The function generating £(s, t) can be computed sim- 
ilar to Eqs. (12 1 and (14 1. In this case, it is a func- 
tion of two variables: h% (x, y) = J2 S t>o £( s > t)y Sxt - Two 
more distributions are needed: S(s, t) and T(s, t) are the 
probabilities of arriving at a vertex of the given degree 
(or degrees) which is at an odd or even distance from 
the starting vertex, respectively, and which belongs to 
a cluster of s extra vertices at odd distance, and t at 
even distance. The corresponding generating functions 
depend on each other: 

hs{z, y) = i - E u i r i + y E n 9 r 9 [M x > y)] q l > 



(Al) 



hr{x, y) = 1 - ^2 n ? r 9 + x ^2 n ? r 9 [hs(x, y)} 



9-1 



(A2) 



and the function generating £(s,t) is 

h( x >v) = I ^ n ? i h s{x,y)] q . 



(A3) 



This allows to compute £(s, t) by taking partial deriva- 
tives in x and y. As in the case of Eqs. (14) and (19 1, 



h^(x, y) is in general a transcendental function and has to 
be solved numerically. However, in some cases it can be 
solved analytically. In the case of 2-swap only (LT2 = 1, 



n g ^2 = 0), Eq. (A3 1 simplifies to the closed form 
x{l-ri) 2 {l + r iy ) 2 



h(.{x,y) 



(1 — r\xy) 2 



(A4) 



The probability £(s, t) in Eq. (22) is then 



1 d'd'h^x^y) 
s\t\ dy s dx t 
2 

1 + s-t 



x,y—Q 

{i-nfrl +t - 1 



t 



(A5) 



if |s — t\ < 1 and otherwise. 

Alternatively, for the case of a single target degree, 
^(rand) ^ s0 computed exactly up to the n-th order 



in r 9 _i by the branching process depicted in Figure 18 
The process begins at step 0, with fco = 1 vertices of 
degree q. At step 1, ki vertices out of qk = q are of 
degree q with binomial probability 



At following steps i > 2 in the branching process, there 
are (q—1) new vertices for each previous vertex of degree 
q. Thus, in every step, k{ vertices are of degree q with 
probability 



f(q- l)fci_A fc 



(g-l)fc < _i-fc i> 



Operations are made on vertices at even steps. Note that 
every new step in the branching process involves higher 
orders in r q -i- Therefore, the expansion of r\ q up to order 
n is obtained by summing the contributions of the first 
n steps: 



x rfj^ k '(l - r,.!)^'- 2 )^ *«)-*», (A6) 

where the sum in {ki} sums for fc = 1, k% — 0, 1, . . . , qk a 
and k t > 2 = 0, 1, . . . , (q - 
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